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Abstract 

We present a computational scheme for ab initio total-energy calculations of materials with 
strongly interacting electrons using a plane- wave basis set. It combines ab initio band structure 
and dynamical mean-field theory and is implemented in terms of plane-wave pseudopotentials. The 
present approach allows us to investigate complex materials with strongly interacting electrons and 
is able to treat atomic displacements, and hence structural transformations, caused by electronic 
correlations. Here it is employed to investigate two prototypical Jahn- Teller materials, KCuFa 
and LaMnOa, in their paramagnetic phases. The computed equilibrium Jahn- Teller distortion 
and antiferro-orbital order agree well with experiment, and the structural optimization performed 
for paramagnetic KCuFa yields the correct lattice constant, equilibrium Jahn- Teller distortion and 
tetragonal compression of the unit cell. Most importantly, the present approach is able to determine 
correlation- induced structural transformations, equilibrium atomic positions and lattice structure 
in both strongly and weakly correlated solids in their paramagnetic phases as well as in phases 
with long-range magnetic order. 

PACS numbers: 71.10.-w, 71.15.Ap, 71.27.-ha 
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I. INTRODUCTION 



The theoretical understanding of complex materials with strongly interacting electrons is 
one of the most challenging areas of current research in condensed matter physics. Exper- 
imental studies of such materials have often revealed rich phase diagrams originating from 
the interplay between electronic and lattice degrees of freedom.- This makes these com- 
pounds particularly interesting in view of possible technological applications. Namely, the 
great sensitivity of many correlated electron materials with respect to changes of external 
parameters such as temperature, pressure, magnetic and/or electric field, doping, etc., can 
be employed to construct materials with useful functionalities.- 

The electronic properties of materials can be computed from first principles by den- 
sity functional theory in the local density approximation (LDA),- the generalized gradient 
approximation (GGA),-'^ or using the so-called LDA+U method.-*^ Applications of these 
approaches accurately describe the phase diagrams of many simple elements and semicon- 
ductors, and of some insulators. Moreover, they often allow to make correct qualitative 
predictions of the magnetic, orbital, and crystal structures of solids, where the equilibrium 
(thermodynamic) structures are determined by simultaneous optimization of the electron 
and lattice systems.->^>2i>i2i>ii>i^ However, these methods usually fail to describe the cor- 
rect electronic and structural properties of electronically correlated paramagnetic materi- 
als. Hence the computation of electronic, magnetic, and structural properties of strongly 
correlated paramagnetic materials remains a great challenge. 

Here the recently developed combination of conventional band structure theory and dy- 
namical mean-field theory,— the so-called LDA+DMFT computational scheme,— has become 
a powerful tool for the investigation of strongly correlated compounds in both their paramag- 
netic and magnetically ordered states. In particular, it provides important insights into the 
spectral and magnetic properties of correlated electron materials,— especially in the vicin- 
ity of a Mott metal-insulator transition as encountered in transition metal oxides.- Up to now 
implementations of the LDA+DMFT approach utilized linearized and higher-order muffin- 
tin orbital [L(N)MTO] techniques'^ and focused on the investigation of electronic correlation 
effects for a given lattice structure. However, the mutual interaction between electrons and 
ions, i.e., the influence of the electrons on the lattice structure, is then completely neglected. 
LDA+DMFT computations of the volume collapse in paramagnetic Cei^ii^ and Pu2^ and of 
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the collapse of the magnetic moment in MnO^ did include the lattice, but only calculated 
the total energy of the correlated material as a function of the unit cell volume.— In the case 
of more subtle structural transformations, e.g., involving the cooperative Jahn- Teller (JT) 
effect,—!^ the L(N)MTO technique is not suitable since it cannot reliably determine atomic 
positions. This is due to the atomic-sphere approximation within the L(N)MTO scheme, 
where a spherical potential inside the atomic sphere is employed. Thereby multipole contri- 
butions to the electrostatic energy due to the distorted charge density distribution around 
the atom are ignored. Instead, the recently proposed implementation of the LDA+DMFT 
approach, which employs plane-wave pseudopotentials^^'^^'^^'^^'^ and thus avoids the atomic 
sphere approximation, does not neglect such contributions. Thereby it becomes possible to 
describe the effect of the distortion on the electrostatic energy.— 

In this paper, we present a detailed formulation of the LDA+DMFT scheme imple- 
mented with plane-wave pseudopotentials reported earlier.— This scheme allows us to 
compute structural transformations (e.g., structural phase stability and structure optimiza- 
tion) caused by electronic correlations. Most importantly, it is able to determine correlation- 
induced structural transformations in both paramagnetic solids and long-range ordered 
solids. Therefore, the present computational scheme overcomes the limitations of stan- 
dard band-structure approaches and opens the way for fully microscopic investigations of 
the structural properties of strongly correlated electron materials. 

We apply this method to study orbital order and the cooperative JT distortion in two 
prototypical JT materials, KCuFs and LaMnOa, and compute the electronic, structural, 
and orbital properties in their room-temperature paramagnetic phase. At room tempera- 
ture, both compounds have a distorted perovskite structure with a strong cooperative JT 
distortion. Considering this structural complexity, KCuFs has a relatively high (tetragonal) 
symmetry,— where only the cooperative JT distortion of the CuFg octahedra is responsible 
for the deviation from the cubic perovskite symmetry. LaMnOa instead crystallizes in a 
more complex (orthorhombic) structure in which MnOe octahedra are simultaneously JT 
distorted and tilted with respect to the ideal cubic perovskite structure.— "^"^ The JT dis- 
tortion persists up to the melting temperature ~ 1000 K in KCUF3. By contrast, in LaMnOa 
it persists only up to Tjt ~ 750 K, the temperature at which the JT distortion vanishes and 
where LaMnOs undergoes a structural phase transition with a volume collapse to a nearly 
cubic structure without JT distortion and orbital order.— 1^ Concerning the electronic con- 
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figuration, KCuFs nominally has a Cu^"*" 3(P structure, i.e., a single hole in the 3d shell. 
By contrast, due to Hund's rule coupling LaMnOs has a single e^f electron bound to the 
fully spin-polarized three t2gt electrons in a high spin 3d* (^291^9?) electronic configuration. 
To properly describe this correlated state, a different treatment compared to KCuFs is re- 
quired, which takes into account the effective on-site spin interaction between and Cg 
electrons.— At low temperatures, both systems display A-type long-range antiferromag- 
netic order, consistent with the Goodenough-Kanamori- Anderson rules for a superexchange 
interaction with antiferro-orbital order. In both compounds the Neel temperature (T/v ~ 38 
K in KCuFs and T/v ~ 140 K in LaMnOs) is remarkably lower than Tjt- 

In this paper we will show that our approach can explain the orbital order, cooperative 
JT distortion, and related structural properties in both materials, in spite of their chemical, 
structural, and electronic differences. The scheme is robust and makes it possible to address, 
on the same footing, electronic, magnetic, and structural properties of strongly correlated 
materials. 

The paper is organized as follows. In Section II we present computational details needed 
to reproduce the results of our calculations. The crystal structures, magnetic properties, 
and results of electronic structure calculations of paramagnetic KCuFs and LaMnOs are 
presented in Sections III and IV, respectively. Finally, the results are summarized in Section 
V. 

II. COMPUTATIONAL DETAILS 

In order to compute the electronic structure of correlated electron materials, we have 
implemented DMFT within a realistic electronic structure approach, which is formulated in 
terms of plane-wave pseudopotentials.— 1^ Following the paper of Anisimov et alM- and 
Trimarchi et al.^ this can be achieved by applying a projection onto atomic-centered 
symmetry-constrained Wannier orbitals,— which gives an effective low-energy Hamiltonian 
Hdft for the partially filled correlated orbitals (e.g., 3d orbitals of transition metal ion). 

The Hamiltonian Hdft provides a realistic description of the material's single-electron 
band structure. It is supplemented by on-site Coulomb interactions for the correlated or- 
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bitals, resulting in a many-body Hamiltonian of the form: 

im 

+ ^ {V - 6^^iJ)nimanim'a' - Hdc- (1) 
im^m' era' 

Here the second and third terms on the right-hand side describe the local Coulomb inter- 
action between electrons in the same and in different correlated orbitals, respectively, with 
V = U — 2J , and Hoc is a double counting correction which accounts for the electronic in- 
teractions already described by DFT (see below). The Coulomb repulsion U and Hund's rule 
exchange J can be evaluated using a constrained DFT scheme within a Wannier-functions 
formalism, making the Hamiltonian (1) free of adjustable parameters.— 

The many-body Hamiltonian (1) is then solved by dynamical mean- field theory (DMFT)^ 
with the effective impurity model treated, e.g., by the numerically exact Hirsch-Fye quantum 
Monte-Carlo method.— Finally, applying a maximum entropy treatment of Monte-Carlo 
data, one obtains the real-frequency spectral functions, which can be further compared to 
physically observable spectra. 

The total energy is another important quantity which can be evaluated from 
DFT+DMFT calculation using the following expression^>2^ 

E = EdfAp] + {Hdft) - J2 'Z^ + i^u) - Enc- (2) 

m,k 

Here Edft[p] is the total energy obtained by DFT. The third term on the right-hand side 
of Eq. (2) is the sum of the DFT valence-state eigenvalues which is evaluated as the thermal 
average of the DFT Hamiltonian with the non-interacting DFT Green's function G^^'^ {iuOn)'- 

E 'Z"" = ^ E Tr[HnFT{^^)G^'\^u^n)]e--'\ (3) 

m,k iiOn,^ 

In this expression, we have assumed that the DFT total energy has only a weak temper- 
ature dependence via the Fermi distribution function, i.e., one neglects the temperature 
dependence of the exchange-correlation potential. {Hoft) is evaluated similarly but with 
the full Green's function including the self-energy. To calculate these two contributions, the 
summation is performed over the Matsubara frequencies iuj^ (typically with n^nax = 10'^), 
taking into account an analytically evaluated asymptotic correction (see below). Thus, for 
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(Hdft) one has 

= T Yl Tr{HnFT{k)[G^{tu^^) - j^]} 

tLU„,K 

+ Ij^'^riHnFTm - -^J2'^r[HnFAk)m^] (4) 

k k 

where the first moment is computed as = HdftO^) + S(ioo) — /i; the asymptotic 
part of the self-energy S(zcxd) is calculated as the average of over the last several iUn 

points. The interaction energy (Hu) is computed from the double occupancy matrix. The 
double-counting correction Edc is evaluated as the average Coulomb repulsion between the 
Nd correlated electrons in the Wannier orbitals. In the case of a paramagnet it corresponds 
to Euc = \UNfi{Nd — 1) — \ JN(i{N(i — 2). Since the Hamiltonian involves only correlated 
orbitals the number of Wannier electrons A'^^ is conserved. Therefore, the double-counting 
correction amounts to an irrelevant shift of the total energy. 

Within this approach, we can determine correlation-induced structural transformations, 
as well as the corresponding change of the atomic coordinates and of the unit cell shape. 
The result can be further used to explain the experimentally observed structural data and 
to predict structural properties of real correlated materials. In Sections III and IV we will 
apply this method to determine the orbital order and the cooperative JT distortion in the 
paramagnetic phase of two prototypical JT systems KCuFa and LaMnOs. 



III. APPLICATION TO KCuFg 

A. Crystal structure and magnetic properties 

KCuFs is the prototype of a material with a cooperative Jahn- Teller (JT) distortion, 
orbital order, and low-dimensional magnetism.— At room temperature it crystallizes in a 
pseudocubic perovskite crystal structure^ which is related to the crystal structure of high- 
Tc superconductors and colossal magnetoresistance manganites and, particularly, to their 
parent compound, LaMnOs. Due to the particular orbital order in KCuFs it is one of the 
rare examples of an ideal one-dimensional antiferromagnetic Heisenberg system.— Thus, the 
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copper ions have an octahedral fluorine surrounding and are nominally in a Cu^"*" {3d^) 
electronic configuration with a single hole in the eg states. This gives rise to a strong JT 
instability that lifts the cubic degeneracy at Cu Cg states due to a cooperative JT distortion.— 
The latter is characterized by CuFq octahedra elongated along the a and b axis and arranged 
in an antiferro-distortive pattern in the ab plane.— At room temperature, there are two 
different structural polytypes with antiferro (a-type) and ferrolike (rf-type) stacking of the 
ab planes along the c axis.— The JT distortion is associated with the particular orbital order 
in KCuFa, in which a single hole alternatingly occupies d^^-z^ and dy2_^2 orbital states, 
resulting in a tetragonal compression (c < a) of the unit cell. The mechanism responsible for 
the orbital order in KCuFs is still being debated in the literature.->2ii2^i^i^i^i^ In particular, 
purely electronic effects such as in the Kugel-Khomskii theory2i and the electron-lattice 
interaction^ have been discussed as possible driving forces behind the orbital order. 

The relatively high (tetragonal) symmetry makes KCuFa one of the simplest system to 
study. In particular, in order to describe the JT distortion, only a single internal structure 
parameter, the shift of the in-plane fluorine atom from the Cu-Cu bond center, is needed. 
Moreover, KCuFs has a single hole in the "id shell resulting in absence of multiplet ef- 
fects. Altogether, the electronic and structural properties of KCuFa have been intensively 
studied by density functional theory in the local density approximation (LDA),- the gener- 
alized gradient approximation (GGA),-"^ or using the so-called LDA+U approach.-i^ While 
the LDA+U calculations account rather well for the value of equilibrium JT distortion in 
KCuFa,— the calculations simultaneously predict long-range antiferromagnetic order which 
indeed occurs in KCuFs below T/v (~ 38 K and 22 K for the a polytype and for the d 
polytype, respectively).— Note, however, that the Neel temperature is much lower than the 
critical temperature for orbital order which is generally assumed to be as large as the melting 
temperature (~ 1000 K). The LDA+U calculations give the correct insulating ground state 
with the long-range A-type antiferromagnetic and d^^-z'^ / dy2_z2 antiferro-orbital order 
consistent with the Goodenough-Kanamori- Anderson rules for a superexchange interaction. 
Non-magnetic LDA calculations instead predict a metallic behavior. Moreover, the elec- 
tronic and structural properties of KCuFa have been recently reexamined by means of 
LDA+U molecular-dynamic simulations, indicating a possible symmetry change and chal- 
lenging the original assignment of tetragonal symmetry.— This symmetry change seems to 
allow for a better understanding of Raman,— electronic paramagnetic resonance,—"^ and 
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x-ray resonant scatterings^ properties at T ^ T^. However, the details of this distortion 
have not been fully resolved yet. 

The LDA+U approach is able to determine electronic properties and the JT distortion in 
KCuFs rather well,— but the application of this approach is limited to temperatures below 
T/v. LDA+U cannot explain the properties at T > T/v and, in particular, at room tempera- 
ture, where KCuFa is a correlated paramagnetic insulator with a robust JT distortion which 
persists up to the melting temperature. 

Here we present an application of the GGA+DMFT computational scheme formulated 
in terms of plane-wave pseudopotentials^^"^ which allows us to determine the structural 
properties, in particular, the amplitude of the equilibrium JT distortion and its energetics, 
in paramagnetic KCuFs. We also report results of a structural optimization, at constant 
volume (constant external pressure) and lattice symmetry, including optimization of the unit 
cell shape and fluorine atomic positions. 

B. Electronic structure and orbital order 

In this section, we present results of the GGA and GGA+DMFT electronic structure 
calculations of paramagnetic insulator KCUF3. In these calculations we have used the ex- 
perimental room-temperature crystal structure with space group I A /mem and lattice con- 
stants a = 5.855 and c = 7.852 A.— The calculations were performed for different values of 
the JT distortion defined accordingly as 5jt = \{di — ds)/{di + ds) where di and ds denote 
the long and short Cu-F bond distances in the ab plane of CuFg octahedra, respectively, 
and 2{di + dg) = a. In the following we express the distortion 6jt in percent of the lattice 
constant a, e.g., 5jt = 0.002 = 0.2 %. In our investigation we consider 0.2% < 6jt < 7%. 
The structural data^^ at room-temperature yield 6jt = 4.4 %. In the present calculations 
we keep the lattice parameters a and c and the space group symmetry fixed, whereas the 
structural optimization involving change of both the JT distortion and lattice constants will 
be discussed in the following section. 

We first calculate the non-magnetic GGA electronic structure of KCuFs, employing the 
plane-wave pseudopotential approach.-i^ For these calculations we use the Perdew-Burke- 
Ernzerhof exchange-correlation functional together with Vanderbilt ultrasoft pseudopoten- 
tials for copper and fluorine, and a soft TrouUier-Martin pseudopotential for potassium. The 
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nonlinear core correction to the exchange-correlation potential has been included to account 
for the overlap between the valence and core electrons. All calculations are carried out in a 
tetragonal unit cell which contains two formula units (10 atoms) per primitive unit cell. We 
use a kinetic-energy cutoff of 75 Ry for the plane-wave expansion of the electronic states. 
The integration in reciprocal space is performed using a [8,8,8] Monkhorst-Pack A;-point grid. 

For all values of (5jt considered here, the non-magnetic GGA yields a metallic rather than 
the experimentally observed insulating behavior, with an appreciable orbital polarization due 
to the crystal field splitting. This is shown in Fig. [1] which presents the GGA density of 
states and the corresponding band structure calculated for 6jt = 4.4 %. Overall, the non- 
magnetic GGA results qualitatively agree with previous band-structure calculations.-'^'^ 
The Cu t2g states are completely occupied and located at about 1-2 eV below the Fermi 
level. Partially filled bands at the Fermi level originate from the Cu eg orbitals. We note 
that an increase of the JT distortion results in a considerable enhancement of the crystal 
field splitting between — y"^ and — bands (in the local framed) that leads to an 
appreciable depopulation of the x"^ — orbital. There is a minor narrowing of the t2g and eg 
bands of ~ 0.2 and 0.1 eV, respectively, as well as a slight up-shift of the center of gravity 
of the t2g bands (~ 0.1 eV) with decreasing JT distortion. 

In Fig. [2] we display our results for the GGA total energy as a function of the JT dis- 
tortion (5jT. Notice that, in agreement with previous studies,-'^ the electron-lattice inter- 
action alone is found insufficient to stabilize the orbitally ordered insulating state. The 
non-magnetic GGA calculations not only give a metallic solution, but its total energy pro- 
file is seen to be almost constant for (5jt < 4 %, with a very shallow minimum at about 
2.5 %. This would imply that KCuFs has no JT distortion for temperature above 100 K, 
which is in clear contradiction to experiment. Obviously, a JT distortion by itself, without 
the inclusion of electronic correlations, cannot explain the experimentally observed orbitally 
ordered insulating state in paramagnetic KCuFs. 

To proceed further, we consider the partially filled Cu eg orbitals as correlated orbitals 
and construct an effective low-energy Hamiltonian Hdpt for each value of the JT distortion 
(5jT considered above. This is achieved by employing the pseudopotential plane-wave GGA 
results and making a projection onto atomic-centered symmetry-constrained Cu eg Wannier 
orbitals.— The resulting Cu x"^ — y"^ and — Wannier orbitals calculated for 5jt = 4.4 % 
are shown Fig. [31 Taking the local Coulomb repulsion U = 7 eV and Hund's rule exchange 
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FIG. 1: (color online) Band structure and orbitally resolved Cu t2g and Cg spectral densities of 
paramagnetic KCuFa as obtained by non-magnetic GGA for 5jt = 4.4%. The zero of energy 
corresponds to the Fermi level. 
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FIG. 2: (color online) Comparison of the total energies of paramagnetic KCuFs computed by GGA 
and GGA-I-DMFT(QMC) as a function of the JT distortion. Error bars indicate the statistical 
error of the DMFT(QMC) calculations. 

J = 0.9 eV~ into account, we obtain the many-body low-energy Hamiltonian (1) for the two 
(m = 1, 2) Cu Cg orbitals, which is further solved (for each value of (5jt) within the single-site 
DMFT using Hirsch-Fey quantum Monte Carlo (QMC) calculations.—"^ The calculations 
have been performed at T = 1160 K (/5 = 10 eV~^), using 40 imaginary-time slices. The 
Matsubara sums in Eq. 2 have been taken over rimax = 10^ frequencies; this gives accuracy 
in the total energy calculation better than 10 meV per formula unit. 

Using now the expression of Eq. [21 we have calculated the GGA+DMFT total energy 
for each value of the JT distortion 5jt considered here. The result of the paramagnetic 
GGA-I-DMFT computation of the total energy is presented in Fig. [2], where it is compared 
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FIG. 3: (color online) lA/mcm primitive cell and Wannier Cg orbitals (x^ — and 3z'^ — in 
left and right, respectively) as obtained by non-magnetic GGA for KCuFs with (5jt = 4.4 %. The 
fluorine atoms and fluorine octahedra are shown in green, the potassium in red, and the Wannier 
Cu Cg charge density in blue. The local coordinate system is chosen with the z direction defined 
along the longest Cu-F bond of the CuFg octahedron. 

with the non-magnetic GGA calculation. In contrast to the GGA result, the inclusion of the 
electronic correlations among the partially filled Cu eg states in the GGA+DMFT approach 
leads to a very substantial lowering of the total energy by ~ 175 meV per formula unit. This 
implies that the strong JT distortion persists up to the melting temperature (> 1000 K), in 
agreement with experiment. This finding is in strong contrast to the absence of JT distortion 
above 100 K predicted by GGA. The minimum of the GGA+DMFT total energy is located 
at the value Sjt ~ 4.2 %, which is also in excellent agreement with the experimental value 
of 4.4 %.— Note however that the total energy minimum position depends on the value of 
Coulomb interaction parameter U. Thus, the calculations of the total energy minima for 
U = 6 eV and U = 8 eV result in optimal JT distortions of 4.15 % and 4.6 %, respectively. 
We note that GGA+DMFT calculations correctly describe both electronic and structural 
properties of paramagnetic KCuFa. This shows that the JT distortion in paramagnetic 
KCuFa is caused by electronic correlations. 

Figure H] shows the spectral density of paramagnetic KCuFs, obtained from the QMC 
data by the maximum entropy method for several values of the JT distortion 6jt Most im- 
portantly, a paramagnetic insulating state with a substantial orbital polarization is obtained 
for all (5jT considered here. The energy gap is in the range 1.5-3.5 eV, and increases with 
increasing of 5jt. The sharp feature in the spectral density at about —3 eV corresponds 
to the fully occupied 3^^ — orbital,— whereas the lower and upper Hubbard bands are 
predominantly of — character and are located at —5.5 eV and 1.8 eV, respectively. The 
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FIG. 4: (color online) Orbitally resolved Cu eg spectral densities of paramagnetic KCuFs as ob- 
tained by GGA+DMFT(QMC) for different values of the JT distortion. The resulting orbitally 
resolved spectral density which is shown here by solid [dashed] line is predominantly of — 
[Sz^ — r^] character (in the local framed). 

corresponding Cu Cg Wannier charge density calculated for the experimental value of JT dis- 
tortion of 4.4 % is presented in Fig.O The GGA+DMFT results clearly show an alternating 
occupation of the Cu d^2_z2 and dy2_z2 hole orbitals, corresponding to the occupation of a 

— hole orbital in the local coordinate system,— which implies antiferro-orbital order. 

The above calculations have been performed for the paramagnetic phase of KCUF3. The 
Neel temperature (T/v ~ 22 — 38 K for different types of orbital order—) is much lower than 
the temperature of present calculations. However, it is known that the ordering temper- 
ature might be overestimated by the single-site DMFT approximation, as is common for 
mean-field theories. To prove the stability of the paramagnetic solution at high tempera- 
tures (with respect to the A-type antiferromagnetic one) we have carried out spin-polarized 
GGA+DMFT calculation at T = 560 K. This calculation has been performed for the A-type 
antiferromagnetic structure using the experimental room-temperature crystal structure of 
KCuFa with 5jt = 4.4 %. However, in agreement with experiment, the calculation gives 
paramagnetic insulating solution with the orbital order as it has been found above. 
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FIG. 5: (color online) I A/ mem primitive cell and hole orbital order as obtained by the 
GGA+DMFT calculation for paramagnetic KCuFa with (5jt = 4.4 %. The fluorine atoms and 
fluorine octahedra are shown in green, the potassium in red, and the Wannier Cu eg charge density 
in blue. The local coordinate system is chosen with the z direction defined along the longest Cu-F 
bond of the CuFg octahedron. 

C. Optimized structure 

To proceed further, we perform a structural optimization of paramagnetic KCUF3. For 
simplicity, the optimization was performed only for two independent structural parameters, 
the lattice constant a and the JT distortion 5jt, keeping the space group symmetry (tetrag- 
onal I A /mem) and the experimental value of the unit cell volume (taken at the ambient 
pressure at room temperature) unchanged.— 

The calculations have been performed in two steps. In the first, we calculate non-magnetic 
GGA electronic structure for different values of (5jt and lattice constant a. Note that in order 
to keep the unit cell volume constant, the value of tetragonal distortion (c/a) was changed 
accordingly. In Figure [6] (a) we plot the total energies obtained by GGA for different 
JT distortion 5jt as a function of the lattice constant a. The data points were further 
interpolated by smooth curvies, whereas the result of the total energy variation — the line 
that connects the minima of the corresponding curves — is marked by black arrows. We 
note that the result of the GGA structural optimization, the variation of the total energy, 
is seen to be constant for (5jt < 2 % with the end point at a ~ 5.75 A. This implies the 
absence of the cooperative JT distortion and results in a nearly cubic (c/a 1.0) unit cell, 
which is in clear contradiction to experiment.— 

In the second step, we construct the effective low-energy Hamiltonian for the partially 
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FIG. 6: (color online) Comparison of the total energies of paramagnetic KCuFs computed by GGA 
(a) and GGA+DMFT(QMC) (b) for different values of the JT distortion 6jic as a function of the 
lattice constant a. Result of the total energy variation is marked by black arrows. Error bars 
indicate the statistical error of the DMFT(QMC) calculations. 

filled Cu Cg orbitals for each value of the JT distortion 5jt and the lattice constant a 
considered here, and compute the corresponding total energies using GGA+DMFT approach 
(see Fig.E](b)).— In contrast to the structural optimization within GGA, the inclusion of the 
electronic correlations among the partially filled Cu eg states in the GGA+DMFT method 
not only correctly describes the spectral properties, but also leads to a very prominent 
minimum in the resulting total energy variation. The minimum is located at the value 
a = 5.842 A and 5jt ~ 4.13 %, which is in excellent agreement with experimental value 
a = 5.855 A and (5jt ~ 4.4 %. Note that in contrast to GGA, the structural optimization 
within GGA+DMFT also correctly predicts the tetragonal compression of the unit cell with 
c/a ^ 0.95.^ 
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IV. APPLICATION TO LaMnOg 



A. Crystal structure and magnetic properties 

LaMnOg is another prototype of a material with a cooperative Jahn- Teller distortion and 
orbital order. It stands in line with the colossal magnetoresistance manganites whose parent 
compound it is.— At ambient pressure and temperature it has an orthorhombic GdFeOg-like 
crystal structure with space group Pnma and four formula units (20 atoms) per primitive 
cell.— The Mn ions have octahedral oxygen surrounding and are in a high-spin 3d^ electronic 
configuration due to Hund's rule coupling, with three electrons in the t2g^ orbitals and a 
single electron in an e^-f orbital (^23^9 orbital configuration). There are two types of structural 
instabilities which give rise to the changes relative to the cubic perovskite structure. The 
first is a JT instability due to the orbital degeneracy that lifts the cubic degeneracy at Mn eg 
states due to developing the cooperative JT distortion of the MnOg octahedra. The second is 
related to a large ion-size misfit parameter \/2{Ro + Rmu) / {Ro + Rlo) which favors rotations 
of the octahedra to accommodate a more efficient unit cell space filling. Rmhi Rbai and Rq 
are the ionic radii of Mn, La, and O ions, respectively. The cooperative JT distortion lifts the 
Cg-orbital degeneracy and leads to an alternating occupation of d^^'^^r'^ and d^y2_j.2 electron 
orbitals in the ah plane (antiferro-orbital ordering) and to a tetragonal compression of the 
unit cell. The rotations of the octahedra lower the symmetry further, finally leading to the 
orthorhombic unit cell. In the paramagnetic Pnma phase the JT distortion experimentally 
persists up to Tjt ^ 750 K. At this temperature, LaMnOg undergoes a structural phase 
transition,— with volume collapse^^ to a nearly cubic structure in which orbital order and 
JT distortion vanish.— A quenching of the JT distortion has also been reported in the room- 
temperature paramagnetic phase under hydrostatic pressure ~ 18 GPa. At ~ 32 GPa it is 
followed by an insulator-metal transition.— 

At temperatures T < T/v ~ 140 K, which are much lower than Tjt ~ 750 K where the 
JT distortion vanishes, LaMnOs shows A-type long-range antiferromagnetic order consis- 
tent with the Goodenough-Kanamori-Anderson rules for a superexchange interaction with 
dj,x'^_T.2j d-iyi_^2 antiferro-orbital order.— 1^ This is also found in spin polarized LDA/GGA 
and LDA-I-U calculations using the experimental values of the crystal structure parameters. 
In this particular case, both magnetic LDA/GGA and LDA+U calculations result in the 
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qualitatively correct insulating ground state with long-range A-tjpe antiferromagnetic and 
antiferro-orbital order.— However, the subsequent structural optimization within the mag- 
netic LDA/GGA calculations results in a metallic solution with reduced JT distortion.— In 
this situation, only the LDA+U scheme is found to give, at equilibrium, the correct insulating 
character of the low-temperature antiferromagnetic phase and a JT distortion in satisfac- 
tory agreement with experiment.— Nevertheless, we have to note again that application of 
this approach is limited to temperatures below T/v. Therefore, LDA+U cannot describe the 
properties of LaMnOa at T > T/v and, in particular, at room temperature, where LaMnOa 
is a correlated paramagnetic insulator with a robust JT distortion. The electronic properties 
of paramagnetic LaMnOa have already been studied within the LDA+DMFT approach. In 
particular, Pruschke and Zolfl^ studied the electronic and magnetic properties and found 
an additional increase of the orbital polarization below T/v. Yamasaki et alM examined 
the electronic structure in order to address the origin of the high-pressure metal-insulator 
transition. Pavarini and Koch^ investigated the temperature dependence of the orbital po- 
larization to find the origin of the cooperative JT distortion and orbital order. However, no 
attempt has been made to determine the structural properties and, in particular, the value 
of the cooperative JT distortion of paramagnetic LaMnOs so far. 

We present here the results of an application of the GGA+DMFT computational scheme 
formulated in terms of plane- wave pseudopotentials^^i^^ to study the electronic and structural 
properties of paramagnetic LaMnOs. These are the first results of a structural optimization 
where the stability of the cooperative JT distortion in paramagnetic LaMnOs was investi- 
gated. In principle, this application can be further extended to investigate the structural 
stability as a function of temperature. Such a full structural optimization will be interesting 
to study the disappearance of the JT distortion at T ~ Tjt- However, this is beyond the 
scope of the present work. 

B. Electronic structure and orbital order 

In this section, we turn to the results of the GGA and GGA+DMFT electronic structure 
calculations of paramagnetic LaMnOs. In these calculations, we have used the orthorhombic 
Pnma crystal structure as reported by Elemans et al, with lattice constants a = 5.742, 
b = 7.668, and c = 5.532 A.— Similar to KCuFs we change the value of JT distortion (5jt, 
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FIG. 7: (color online) Orbitally resolved Mn t2g and Cg spectral densities of paramagnetic LaMnOs 
as obtained by non-magnetic GGA for 5jt = 0.138. The zero of energy corresponds to the Fermi 
level. 

which is now defined as the ratio between the difference of the long (di) and the short {dg) 
bond distances and the mean Mn-0 distance in the basal ab plane, i.e., 5jt = 2{di — ds) / {di + 
dg). Structural datai^i yield 5jt = 0.138. Note that in the calculation we change only the 
parameter 5jt (0 < 5jt < 0.2) and keep the value of the MnOe octahedron tilting and 
rotation fixed. 

Using the plane-wave pseudopotential calculation scheme, we calculated the non-magnetic 
GGA electronic structure^"^ of paramagnetic LaMnOs. We employed the Perdew-Burke- 
Ernzerhof exchange- correlation functional together with Vanderbilt ultrasoft pseudopoten- 
tials, including a nonlinear core correction to the exchange-correlation potential. All calcula- 
tions were carried out in a 20-atom orthorhombic Pnma unit cell. We used a kinetic-energy 
cutoff of 45 Ry for the plane-wave expansion of the electronic states. The integration in 
reciprocal space was performed using a [10,10,10] Monkhorst-Pack /c-point grid. 

For all values of 5jt considered here the non-magnetic GGA calculations give a metallic 
solution with a considerable orbital polarization due to the crystal field splitting. Overall, 
these results qualitatively agree with previous band structure calculations,— namely, that 
the GGA cannot describe a paramagnetic insulating behavior which is found in experiment. 
We notice that even for the large 5jt value of 0.2 (~ 45 % larger than found in experiment^) 
the GGA calculations predict a metal. The non-magnetic GGA density of states and the 
corresponding band structure calculated for the JT distortion 6jt = 0.138 are presented in 
Fig. [71 In contrast to KCuFs, the partially filled bands at the Fermi level now originate 
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FIG. 8: (color online) Comparison of the total energies of paramagnetic LaMnOs computed by GGA 
and GGA+DMFT(QMC) as a function of the JT distortion. Error bars indicate the statistical 
error of the DMFT(QMC) calculations. 

from the Mn t^g orbitals. This is due to the Sc?^ electronic configuration of Mn ions and 
the fictitious paramagnetic state without local moments obtained by non-magnetic GGA. 
In fact, the Hund's rule coupling results in a strong on-site spin polarization of the Mn t2g 
and Cg orbitals. Therefore, the t2g\ states become completely occupied and located below 
the Fermi level, while the remaining electron fills the tg^ states. An increase of the JT 
distortion results in considerable enhancement of the crystal field splitting between — y'^ 
and 32;^ — bands (in the local framed) which reaches ~ 1.1 eV for 5jt = 0.138. The overall 
Cg band width is about 2.8-3.0 eV, which is remarkably much smaller than the estimates of 
the Coulomb interaction parameter U found in the literature.—'^ 

In Fig. [8] we display our results for the GGA total energy as a function of the JT dis- 
tortion 5jT. In contrast to experiment, the non-magnetic GGA calculations give a metallic 
solution without cooperative JT distortion. Thus, the GGA total energy is almost parabolic 
which implies the absence of a cooperative JT distortion and is in clear contradiction to 
experiment.— 1^ As in the case of KCuFs this shows the importance of electronic corre- 
lations, without which the experimentally observed orbitally ordered, insulating state in 
paramagnetic LaMnOs cannot be explained. 

Next, we turn to the GGA-I-DMFT results where we treat the Mn Cg orbitals as correlated 
orbitals. Using pseudopotential plane-wave approach, we perform a projection onto atomic- 
centered symmetry-constrained Mn Cg Wannier orbitals,— which are shown in Fig. [9l In 
this calculation we assume that three (among the 3(i^ electronic configuration) electrons are 
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FIG. 9: (color online) Pnma primitive cell and the Wannier Cg orbitals {x'^ — and — in 
left and right, respectivel}«^) for LaMnOs with 6jt = 0.138 according to the non-magnetic GGA 
calculation. The oxygen atoms and oxygen octahedra are shown in green, the lanthanum in red, 
and the Wannier Mn Cg charge density in blue. 

localized in the t2g orbitals. Therefore, they are treated as classical spins S, with a random 
orientation above T/v (i.e., there is no correlation between different Mn sites), which couple 
to the Cg electron with an energy JS. This coupling can be estimated as the energy of the 
splitting of the e^j and Cgi bands in the ferromagnetic band-structure calculations and gives 
an additional term in the Hamiltonian (1), namely. 



This corresponds to the ferromagnetic Kondo-lattice model (KLM) Hamiltonian with an on- 
site Coulomb repulsion between Cg electrons, which has been intensively studied as a possible 
microscopic model to explain colossal magnetoresistance in manganites.— We note that 
in order to calculate the total energy one needs to modify Eq. 2 by adding the expectation 
value of the JS term which describes the total energy gain due to the spin polarization of 
the Cg orbitals at the Mn site. 

We take the local Coulomb repulsion U = 5 eV, the Hund's rule exchange J = 0.75 eV, 
and 275" = 2.7 eV from the literature^ and further solve the many-body Hamiltonian (5) 
for each value of 6jt using the single-site DMFT with Hirsch-Fey quantum Monte Carlo 
(QMC) calculations.—"^ The calculations were again performed at T = 1160 K (/5 = 10 
eV~^), using 40 imaginary-time slices. 

In Fig. [8] we present the result of the paramagnetic GGA+DMFT computation of the 



H = Hdft + U 
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FIG. 10: (color online) Orbitally resolved Mn Cg spectral densities of paramagnetic LaMnOs as 
obtained by GGA+DMFT(QMC) for different values of the JT distortion. The resulting orbitally 
resolved spectral density shown by solid [dashed] line is predominantly of 3z^ — [x^ — y^] character 
(in the local framed). 

total energy, where it is compared with the resuhs of the non-magnetic GGA calculation. 
In contrast to the GGA result, the correlated electron problem solved by GGA+DMFT 
approach gives a substantial total energy gain of ~ 150 meV per formula unit. This implies 
that the cooperative JT distortion indeed persists up to high temperatures (T > 1000 K), 
while in GGA a JT distortion does not occur at all. Taking into account that the calculations 
have been performed for the low-temperature crystal structure of LaMnOg^ this estimate 
(150 meV) is in good agreement with Tjt ~ 750 K at which the JT distortion vanishes.—*^ 
However, the structural change as a function of temperature in LaMnOa, as well as the 
disappearance of the orbital-order and JT distortion^"^ remains an open problem. The 
minimum of the GGA+DMFT total energy is located at the value (5jt ~ 0.11, which is also 
in good agreement with the experimental value of 0.138.— We note that GGA+DMFT 
calculations correctly describe both electronic and structural properties of paramagnetic 
LaMnOs. This shows that the JT distortion in paramagnetic LaMnOs is caused by electronic 
correlations. 

The spectral densities of paramagnetic LaMnOs calculated for several values of the JT 
distortion (5jt using the maximum entropy analysis of the QMC data are shown in Fig. [TOl 
For large 5jt, we find a strong orbital polarization which gradually decreases for decreasing 
JT distortion. The occupied part of the Cg density is located at about -2 - -1 eV and corre- 
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FIG. 11: (color online) Pnma primitive cell and orbital order as obtained by the GGA+DMFT 
calculation for paramagnetic LaMnOs with Jjx = 0.138. The oxygen atoms and oxygen octahedra 
are shown in green, the lanthanmn in red, and the Wannier Mn tg charge density in blue. 

sponds to the eg states with spin parallel to the t2g spin at that site. It has predominantly 
Mn d^^2_^.2 and d^y2_^.2 orbital character with a considerable admixture of (1^2 _^2 for small 
JT distortions. The energy gap is about 2 eV for large 5jt and considerably decreases with 
decreasing 5jt, resulting in a pseudogap behavior at the Fermi level for 5jt = 0. In Fig. [11] 
we show the corresponding Mn eg Wannier charge density computed for the experimental 
JT distortion value of (5jt = 0.138. The result clearly show an alternating occupation of the 
Mn d^^2_j.2 and d^y2_r2 orbitals, corresponding to the occupation of a "iz^ — orbital in the 
local frame,— which implies antiferro-orbital order. Thus, in agreement with experiment, the 
calculations give a paramagnetic insulating solution with antiferro-orbital order and stable 
JT distortion. 



V. SUMMARY AND CONCLUSIONS 



In conclusion, by formulating GGA+DMFT — the combination of the ab initio band 
structure calculation technique GGA with the dynamical mean-field theory — in terms of 
plane-wave pseudopotentials,— we constructed a robust computational scheme for the in- 
vestigation of complex materials with strong electronic interactions. Most importantly, the 
computational scheme presented here allows us to explain correlation-induced structural 
transformations, shifts of equilibrium atomic positions and changes in the lattice structure, 
and to perform a structural optimization of paramagnetic solids. We presented apphca- 
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tions of this approach to two prototypical Jahn- Teller materials, KCuFa and LaMnOs, and 
computed the orbital order and cooperative JT distortion in these compounds. In par- 
ticular, our results obtained for the paramagnetic phase of KCuFa and LaMnOa, namely 
an equilibrium Jahn- Teller distortion (5jt of 4.2 % and 0.11, respectively, and antiferro- 
orbital order, agree well with experiment. The present approach overcomes the limitations 
of the LDA+U method and is able to determine correlation-induced structural transforma- 
tions in both paramagnetic and long-range magnetically ordered solids, and can thus be 
employed for the lattice optimization and molecular-dynamic simulations of these systems. 
The GGA+DMFT scheme presented in this paper opens the way for fully microscopic inves- 
tigations of the structural properties of strongly correlated electron materials such as lattice 
instabilities observed at correlation-induced metal-insulator transitions. 
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